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Abstract 


We  studied  the  capability  of  a  Fourier  spectrometer  system  to  remotely  detect  trace  gases 
in  localized  clouds,  e.g.,  stationary  source  effluents.  Detection  is  based  on  the  degree  to 
which  the  observed  IR  spectral  radiance  contrast  between  the  cloud  and  adjacent  background 
is  correlated  with  a  computed  reference  spectrum.  It  is  shown  that  trace  gases  can  be 
reliably  detected  when  spectral  feature.^  are  well  below  the  noise  level.  The  minimum 
detectzLble  quantities  (MDQ's)  for  various  trace  gases  at  one  atmosphere  total  pressure  are 
given.  The  MDQ's  determine  the  combinations  of  gas  column  thickness  and  gas-background 
temperature  difference  that  correspond  to  95  percent  detection  probability  and  one  percent 
false  detection  probability  when  an  FTS  system  with  modest-size  foreoptics  views  the  target 
through  a  path  equivalent  to  approximately  one  air  mass. 


Introduction  ^ 


Infrared-active  trace  gases  in  a  localized  cloud  can  be  remotely  detected  on  the  basis 
of  the  contrast  they  produce  in  an  IR  scene  image.  Of  course,  the  gas  temperature  and 
average  background  brightness  temperature  must  be  different,  and  the  resultant  spatial 
radiance  contrast  must  be  distinguishable  from  background  clutter.  Resolution  of  the  .image 
into  many  spectral  elements  will  give  additional  information  that  can  be  used  to  suppress 
the  effects  of  background  clutter,  atmospheric  spatial  variations  and  sy.'tem  noise,  and  to 
distinguish  between  different  trace  species  in  the  cloud. 

The  imaging  spectrometer  gives  the  spectrum  of  the  difference  between  target  (trace  gas) 
and  background  radiances,  which  can  be  compared  to  a  computed  reference  contrast  spectrum. 
Parameters  based  on  the  degree  of  correlation  between  the  two  spectra  can  be  used  to  decide 
whether  a  particular  target  species  is  present  and  also  to  characterize  its  amount. 

From  studies  reported  elsewhere^  >2  it  was  decided  that  a  mosaic  FTS  was  the  best  type  of 
sensor.  This  choice  was  based  on  an  assumed  target  cloud  size,  detection  range  and  sensi¬ 
tivity  requirements  (foreoptics  size  and  field  of  view),  the  number  and  spectral  locations 
of  the  target  species  IR  bands,  and  the  current  technology  of  high-performance  detector 
mosaics.  This  paper  gives,  for  a  particular  baseline  FTS  system,  the  minimum  detectable 
quantities  of  13  molecular  species.  System  noise-equivalent  spectral  radiances  for  the 
detection  bands  are  also  given,  which  allows  scaling  of  the  MDQ’s  with  respect  to  system 
radiometric  parameters.  The  reported  MDQ's  correspond  to  ideal  conditions;  i.e.,  it  is 
assumed  that  system  performance  is  limited  by  detector  (system)  noise  rather  than 
aurospheric/background  variations,  and  that  the  reference  spectrum  (which  includes  atmos¬ 
pheric  spectral  absorption  effects)  is  known  precisely.  The  effects  of  non-ideal  conditions 
are  discussed  in  Ref.  1. 


Theoretical  basis 

It  can  be  shown^  that  the  measured  spectral  radiance  contrast,  apart  from  system  spectral 
noise,  is 

=  Nj  -Ni  =  DT^Ug^;  D  =  (Tg.T^f^)  (D 

where  "v  '  atmospheric  spectral  transmittance  between  the  target  and  sensor, 

Ug^  '  spectral  absorption  coefficient  of  the  (single)  target  gas, 

u  «  molecular  column  thickness  of  the  target  gas  (number  of  molecules  per  unit 
area  in  the  line-of-sight) ,  and 
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-  differanc*  batwaan  tba  apactral  radlanca  of  a  blackbody  at  tba  gaa  tamparatura 
Tg  and  tha  actual  apactral  radlanca  of  tha  background  at  tamparatura  Tgf ^ , 
avaragad  ovar  tba  dataction  apactral  bond. 

Tha  datactabla  quantity,  danotad  by  0,  la  asaantlally  a  acala  factor  In  the  maaaurad  con¬ 
trast  spectrum.  Ec[uation  (1)  la  an  approximation  only  bacauaa  wa  have  replaced  the  actual 
spectral  variation  ABy  by  its  average  value.  However,  the  variation  is  flow  compared  to 
that  of  Ty  and  Ugy  and  is  nearly  linear  ovar  narrow  dataction  bands,  so  that  tha  approxima¬ 
tion  is  a  good  one,  especially  if  ly  ^  1  and  Sgy  is  mora-or-lass  symmetrical  about  tha 
center  of  the  detection  band. 

A  reasonably  accurate  theoretical  estimate  of  TyOgy  can  be  obtained  using  a  lina-by-line 
computer  code  such  as  AFGL's  FASC0D1.3  Then  an  eatimate  D'  of  D  can  he  obtained  by  finding 
the  value  that  minimizes  the  mean  square  difference  between  the  measured  spectrum  ANy  and 

course,  Ogy  is  a  function  of  the  gas  temperature,  but  31r.Ony/3Tg  <<  31nABy/3Tg. 
Hence,  the  use  of  some  guess  temperature  in  computing  Ugy  will  not  re'jult  in'^a  large  error’ 
in  D' . 


Some  simple  types  of  background  and  atmospheric  interferences  can  be  suppressed  by  sub¬ 
tracting  the  means  of  ANy  and  TyOgy  from  these  spectra  before  they  are  used  to  estimate  D. 
The  least-squares  estimator  of  D  is  then 


]^{a-a)  (b-S) 


(2) 


where  a  i  iNy  +  NOlSEy ,  b  5  Tya_y,  and  the  summations  are  over  the  digitized  spectrum  values. 
Equation  (2)  is  the  prescription  for  computing  the  estimate  D'  gi\en  the  measurement  a  and 
refeyey'.ce  ccrtraat  spectrum  b.  By  combining  Eqs .  (1)  and  (2),  we  find  that  D'  can  also  be 

expressed  by 


D  «  y  - - j 

^b 

where  oj-n  is  the  covariance  of  the  digitized  reference  spectrum  b  and  the  spectrum  noise 
NOISEy,  and  is  ^he  variance  of  the  digitized  reference  spectrum. 

The  spectrum  noise  samples  will  be  Gaussian  with  zero  mean,  and  the  samples  will  also  be 
independent  provided  the  sample  spacing  is  approximately  equal  to  the  resolution  Av  of  the 
FTS  syste.-n.  It  can  be  shown  that  if  the  spacing  is  Av  the  second  term  of  Eq.  (3)  is 


®bn  ,7  NESR 

-  - - -  X - =  XOn-  (4) 

.M  3b  .■M  Ob 

where  x  is  a  unit  normal  random  variate ,  M  is  the  total  num))er  of  resolved  spectral  ele¬ 
ments  in  the  spectra  a  and  b,  and  NESR  is  the  average  syste'n  noise-equivalent  spectral 
radiance  over  the  detection  band.  Note  that  the  effective  rms  noise  ci^  equals  NiSSR 
because  iNy  is  the  difference  of  two  measurements  with  independent  noise.  Thus,  thiH  rms 
uncertainty  in  D'  is  o^-  =  »'T  NESR/(/M  Ob)  • 

The  quantity 


ObD 

rT.  NESR 


.5) 


is  the  rms  signal  spectral  variation  over  the  rms  spectral  noise,  or  the  spectral  signal-to- 
noise.  Since  D/oy.  is  this  quantity  times  yH,  it  is  evidert  that  for  sufficiently  large  M 
the  detectable  quantity  can  be  estimated  accurately  even  ii.  local  spectrum  features  are  wel’ 
below  the  noise  level  (S/N  <<  1)  . 


D  will  be  negative  if  the  target  gas  is  observed  in  absorption  (is  "cooler"  than  the 
background)  and  positive  if  it  is  observed  in  emis-sion.  Hence,  detection  of  the  ta:;:get  gas 
can  be  based  on  the  test 


where  Ry  is  some  predetermined  tlireahold.  Knowing  the  statistics  of  D'  (Eqs.  3  and  4)  we 
can  work  out  the  statistics  for  lD'|  for  both  target  present  (D  f  0)  and  target  not  present 
CD  <•  0)  ,  and  finally  determine  a  threshold  Rp  that  yieltis  satisfactory  detection  and  false 
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detection  probabilltle*.  It  can  be  ahown  that  the  tbie^old  value 

Rd  -  1.58  On-  ■  2.58  (7) 

^  Ob 

resulta  In  approximately  95  percent  detection  probability  and  approximately  one  percent 
false  detection  probability  wh««n  D  equals- 

4.23  o  -  -  4.23  =  mdq  .  (8) 

°  /H  Ob 

This,  by  definition,  Is  the  minimum  detectable  quantity. 

The  NESR,  M  and  o^  all  vary  with  t'.ie  speetroaator  rusolution  Av.  A  near-optimum  value 
for  iv  is  the  halfwidth  of  the  target  gas  spectral  linen,  which  is  -0.1  cm"^  for  targets 
near  sea-level.  That  is,  4v  ~  0.1  cm"l  will  give  the  ssialJest  MDQ  from  Eq.  (8).  The  single 
parameter  in  Eq.  (8)  that  depends  on  the  target  gas  and  detection  geometry  is  o^,  the 
standard  deviation  of  the  reference  spectrum  r degraded  to  0.1  cm'^  resolution. 

FTS  baseline  system 

The  FTS  baseline  configuration  that  was  selected  uses  i  four-port  Michelson  inter¬ 
ferometer  with  a  single  dichroic  2md  two  detector  motiaica  in  each  of  the  two  output  ports 
(see  Fig.  1).  Thir  arrangement  allows  the  simulatneous  .'atection  of  up  to  four  different 
target  species.  The  fourth  port  {a  second  input  port)  provides  for  illumination  of  the 
bac)tside  of  the  interferometer  beamsplitter  by  an  internal  blackbody  reference  or  by  a 
different  (bac)cground)  portion  of  the  external  scene  than  is  viewed  by  the  primary  input 
port.  Utilization  of  the  fourth  port  in  thin  manner  effects  an  optical  subtraction  of  the 
scene  and  reference  blackhody  spectra  or  of  the  displace':'',  scene  spectra,  and,  more  important, 
results  in  a  potentially  lar-ge  reduction  in  the  interfe,-,ogram  dynamic  range.'*  If  the 
reference  blac)cbody  and  background  have  the  same  brightness  temp>erature ,  the  resultant 
target-background  spectral  contrast  obtained  by  this  optical  scheme  is  the  same  as  would  be 
obtained  using  a  conventional  single-port  interferometer  and  merely  subtracting  the  outputs 
of  two  different  mosaic  elements.  The  details  of  the  selected  baseline  configuration  are 
given  in  Ref.  1. 


MDQ  and  fiald  ft  r«Bult» 
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A 


The  MDQ'a  were  calculated  for  several  candidate  detection  bande  of  each  of  14  different 
molecular  species ,  The  "best"  band  for  each  species  was  selected  on  the  ba»’xa  of  the  corres¬ 
ponding  minimum  detectable  column  thicknesses  (MIN  u's)  for  a  given  small  tauperature 
difference  between  target  cloud  and  background.  The  relationship  betmen  u  and  MDU  is, 
of  course,  strongly  dependent  on  the  spectral  location  of  the  detection  ba.id  through  the 
Planck  function  temperature  derivative.  Table  1  is  a  partial  list  of  the  best  bands,  the 
NE^R's  for  a  lO-second  spectrum  scan  time,' the  corresponding  MDQ  and  the  MIN  v.  The  last 
quantity  is  the  minimum  detectable  concentration  in  ppmV  for  the  conditions  noted  in  the 
table  footnote. 


Gas 


HCl 

NOj 

DF 

HBr 

N^O 

CO 

C«4 

SOj 

NH3 

HN03 

^2«4°2 

CH3I 

NO., 


_ Ttible  1.  MDQ  and  MIN  v  reaults _ 

_.  *HeSR  *MDQ  (molec/cm^) 

Band  (cm"'*') _ (W/gm^-sr-cm~^) _ (W/cm^-sr-cm~l) 


2900-3000 

1.3 

X 

10“® 

2.5 

X 

10® 

2650-2935 

1.4 

X 

10'® 

1.3 

X 

0 

rH 

0 

2700-2900 

2.6 

X 

10-5 

3.5 

X 

lo5 

2450-2650 

3.9 

X 

10~5 

2.4 

X 

10^8 

2160-2210 

3.1 

X 

10'5 

3.1 

X 

0 

fH 

0 

2130-2185 

3.6 

X 

10-5 

7.2 

X 

10^° 

1295-1310 

3.9 

X 

10'5 

2.1 

X 

10^2 

1090-1210 

1.5 

X 

10-5 

2.7 

X 

10^1 

915-970 

1.1 

X 

10-8 

2.0 

X 

0 

0 

887.5-902.5 

6 . 0 

X 

10-5 

3.4 

X 

10^“ 

880-1000 

2.0 

X 

10*5 

5. 

X 

10^° 

790-990 

2.4 

X 

10-8 

6.7 

:C 

10^1 

710-795 

1.3 

X 

10-8 

4.0 

X 

10^^ 

**MIN  V  (ppmV) 
20 
92 

J.1 

51 

22 

16 

163 

16 

1.0 

1.6 

2.8 

32 

19 


Note  that  the  MDQ  scales  as  the  system  NESt  (Eq.  8)  .  Tht  more  fundamental  quantity 
for  each  detection  band  can  be  determined  from  Eq.  (8)  and  the  third  and  fourth  columns  of 
Table  1.  The  number  M  is  equal  to  the  width  of  the  band  divided  by  0.1  cm”!.  The  Ojj’s 
and  MDQ's  correspond  to  a  vertical  path  through  most  of  the  Midlatitude  Summer  model 
atmosphero . 

Field  measurements  were  perfomved  at  Utah  State  University  to  val  ate  the  theoreiical 
predictions.  In  the  experiment  a  model  stack  emitting  a  controlled  now  of  N2O  in  front  of 
a  heated  "background"  plate  was  observed  by  an  FTS  system  at  a  distance  of  273  meters.  The 
detection  band  selected  for  t.he  experiment  is  2170-2260  cm" ^ ,  which  includes  most  of  the  V3 
band  of  N2O  (this  is  a  better  N2O  detection  band  than  the  one  given  in  Table  1  for  the 
measurement  geometry  of  the  experiment) .  The  N2O  was  successfully  identified  and  quantified 
in  repeated  trials  when  th«  detectable  quantity  was  slightly  greater  than  the  calculated 
MDQ.  Figure  2  shows  a  computed  (FASCODl)  reference  contrast  spectrum  for  one  of  the  data 
runs,  and  two  examples  of  measured  contrast  spectra  that  resulted  in  N2O  detection.  The 
middle  panel  of  the  figure  shows  the  measured  contrast  when  the  detectable  amount  is  tpproxi- 
m.attly  four  times  the  MOQ;  :.n  the  bottom  panel  the  detectable  amount  is  ~1.5  t  mes  the  MDQ, 
There  is  extremely  high  probability  of  detection  (>>  0,95)  when  the  detectable  amount  is  1.5 
times  the  MDQ,  even  though  the  measured  contrast  spectrum  for  this  amount  looks  like  noise. 
The  spectrum  signal-to-noise  defined  by  Eq.  (5)  is  0.26  for  the  bottom  panel  in  Fig.  2;  at 
the  MDQ  level  the  S/N  would  be  0.15. 
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^For "baseline  system  and  dwell  time  -  10  secs. 

**MIN  V  »  minimum  detectable  volume  concentration  for  t(i  ■ 
10  m,  iT  5'C.  ,  Tgfc  *  300K. 
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Figure  2.  Normalized  reference  contrast  spectrum  for  V3  band  of  N2O  (top  panel) 
and  r.\'o  normalized  contrast  spectra  obtained  from  measurements. 
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